00 

o 

O 
(N 
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It is theoretically shown that the excess liquid-liquid interlacial tension between two electrolyte 
solutions as a lunction ol the ionic strength I behaves asymptotically as 0( — vJ) for small I and 
as C(±7) for large I. The former regime is dominated by the electrostatic potential due to an 
unequal partitioning of ions between the two liquids whereas the latter regime is related to a fi- 
nite interfacial thickness. The crossover between the two asymptotic regimes depends sensitively 
on material parameters suggesting that, depending on the actual system under investigation, the 
experimentally accessible range of ionic strengths can correspond to either the small or the large 
ionic strength regime. In the limiting case of a liquid-gas surface where ion partitioning is absent, 
the image chage interaction can dominate the surface tension for small ionic strength I such that 
an Onsager-Samaras limiting law 0{— Jln(J)) is expected. The proposed picture is consistent with 
more elaborate models and published measurements. 

PACS numbers: 68.05.-n, 82.45. Gj, 68.03. Cd 
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The temporal stability of liquid- liquid emulsions, 
which is of enormous importance for applications in, e. g., 
chemical, pharmaceutical, food, and cosmetic industries, 
largely hinges on the liquid-liquid interfacial tension [l| 
modified by surfactants, cosurfactants, and even colloidal 
particles |2|]. In order to theoretically understand and 
predict the liquid-liquid interfacial tension as a function 
of additives a first step is modeling a liquid-liquid inter- 
face in the presence of electrolytes but in the absence of 
surfactants. Remarkably, the dependence of the liquid- 
liquid interfacial tension on the electrolyte concentration 
is, in contrast to the liquid-gas surface tension [3J], not 
well understood. This is quite astonishing because liquid- 
liquid interfaces have been investigated for a long time by 
means of electrocapillary measurements [J] . The few re- 
ported measurements of the liquid-liquid interfacial ten- 
sion as a function of the ionic strength known to the au- 
thors, Refs. [a, ||j|> [Zj, seem to confirm the linear relation 
at large ionic strengths well-known from liquid-gas sur- 
face tension measurements |8j . At low ionic strengths the 
liquid-gas surface tension exhibits the Jones-Ray effect, 
i. e., a minimum of the surface tension as a function of the 
ionic strength 9] , whose analog for liquid- liquid interfa- 
cial tensions has been addressed in the experimental lit- 
erature, to the authors' knowledge, only in Ref. |£g. The- 
oretical approaches to liquid-gas surfaces are very often 
based on the assumption that the gas phase is completely 
free of ions [3|, which leads to a charge neutral liquid 
phase. Considering the image charge interaction as domi- 
nating the liquid-gas surface tension at low ionic strength 
the Onsager-Samaras limiting law can be derived [101 ] . 
However, assuming a non- vanishing ionic strength in the 
gas phase, Nichols and Pratt found indications that the 
liquid-gas surface tension in some instances can also scale 
with the square root of the ionic strength in the low salt 
limit [ll|. By means of an elaborate Ginzburg-Landau- 
like model for liquid-liquid interfaces, taking ion densities 



and solvent composition explicitely into account, Onuki 
recently observed such a square root behavior for the 
liquid- liquid interfacial tension, too [12J . It is the aim of 
this letter to argue in terms of a minimal model that un- 
equal ion partitioning and charge separation are the key 
features of liquid-liquid interfaces of electrolyte solutions 
at low ionic strength. Onsager-Samaras-like behavior can 
be found only in the absence of unequal ion partitioning 
and is therefore unexpected for liquid-liquid interfaces. 

In order to define the present model consider an infinite 
system composed of two homogeneous solvents A and B 
located within the half spaces z < and z > 0, respec- 
tively, of a Cartesian coordinate system. In the interior 
of the solvents the relative dielectric constant e(z) at po- 
sition z is given by e(z < 0) = ea and e(z > 0) = Sb- 
In the following the abbreviation n :— Wsa/sb will be 
useful. Monovalent ions are distributed in both solvents 
giving rise to local equilibrium number densities g a (z) 
at position z with a = + and a = — denoting cations 
and anions, respectively. Deep in the solvent phases lo- 
cal charge neutrality holds, i. e., g a (—oo) —: qa and 
Q a {po) =: qb- The 'partition coefficient is defined by 
P '■= V ' QAl ' Qb- In general, the solubility of a ions dif- 
fers in the two solvents. This effect can be described 
by solvent-induced potentials V a (z) which take the lim- 
iting values V a (— oo) := and K*(oo) := f a where f a 
is the solvation free energy difference of an a ion in sol- 
vent B as compared to solvent A. Verwey and Niessen 
[13J | assumed the steplike form V^f N (z) — f a Q(z), where 
9 denotes the Heaviside function. Such a model ignores 
interfacial effects due to an actually smooth dielectric 
function e, finite ion size, van der Waals forces, solvation 
(structure making and structure breaking), and image 
charges [3(- All these effects depend on material param- 
eters of the system but they depend, with the exception 
of the image charge interaction, not directly on the ionic 
strength. Moreover, the image charge interaction decays 



as 0(exTp(—2nA,B\z\)/\z\) with kT^b denoting the Debye 
screening length in phase A for z — » — oo and in phase B 

for z — > oo 
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lOL 1 121 ] , whereas the electrostatic potential 
is expected to decay much slower as 0(exp(— «a,b|^|))- 
Hence the image charge interaction is expected to be neg- 
ligible outside the interfacial region. A simple account of 
the mentioned interfacial effects is given by the shifted 
Verwey-Niessen potentials V a (z) :— f a Q(z — s) where 
the discontinuity is located at position z = s, similar 
to the interface model by Johansson and Eriksson [l4j. 
Note that the electrostatic potential is the only interac- 
tion which is not described by the solvent-induced poten- 
tials V a because it is the longest-ranged ionic-strength- 
dependent interaction. Moreover, the shift of the ion 
densities with respect to the solvent composition profile 
in Onuki's work [12| are compatible with the introduction 
of external fields similar to the present solvent-induced 
potentials V a . The location of the discontinuity of the 
solvent-induced potentials with respect to the dielectric 
interface at z = is a property of the solvents and the 
electrolyte. The analysis will in fact reveal that only 
changing the anion type can shift the discontinuity of V a 
to the opposite side of the interface. Without restriction 
s > is assumed, i. e., solvent B is defined as the one 
where the discontinuity of V a is located. 

The equilibrium structure represented by the density 
profiles g a is most easily calculated in terms of density 
functional theory [15| . In units of the thermal energy 
ksT, the elementary charge e, and the vacuum Bjerrum 



length I = 



~knT 



with the permeability of the vac- 



47T£ V ac^ 

uum e V ac, and within a mean-field theory ignoring ion-ion 
correlations, the density functional of the grand potential 
per unit surface area 

ft[Q±] = X! / & z Ba{z)[ \n(g a (zj) - 1 - fi a 



+ V a (z) + a-4>(z,[Q±})) (1) 



is to be minimized with respect to g a . Here \i a is the 
chemical potential of species a and the electrostatic po- 
tential (j)(z, [g±]) at position z, which is a functional of 
the ion density profiles g±, fulfills the Poisson equation 

^e(z)^<Xz, [g±]) = -Me+(z) - Q-(z)) (2) 

with the Dirichlet boundary conditions </>(— oo) = and 
4>{oo) = <f>Di where tfirj '■= |(/- — /+) is the Donnan 
potential following from the local charge neutrality in the 
bulk liquids. The electrostatic potential <j) is continuous 
and it holds ea<J>'{Q~) — es0'(O + ), where a prime denotes 
a spatial derivative. From the Euler-Lagrange equations 
corresponding to Eqs. ([T]) and @ one readily derives for 
the shifted electrostatic potential iJj(z) :— (f)(z) — 4>D^{z — 



s) the linearized Poisson-Boltzmann equation 
d 2 



dz 2 



ip(z) — k(z) ^p(z), z^0,s 



(3) 



with the homogeneous Dirichlet boundary conditions 
■0(±oo) = and the piecewise constant Debye screen- 



ing factor k(z) defined by n(z) 
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for z < 0, k(z) z 
k{z) 2 ■- k% :-- 



k~ := 8irgA/sB for z e (0, s), and 
Sngs/sB for z > s. Moreover, the par- 



tition coefficient is found as p — exp((/ + + /_)/4). The 
solution of Eq. ([3]) is 



D 



exp(K^z) 



,z <0 



ip(z) = < 



_ — — (cosh(Kjz) + nsinh(Kiz)) ,ze(0,s) 
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xp(ncosh(Kis) + sinh(Kis)) ,z>s 



with D := (1 + np) coshes) + (n + p) smh(nis). For the 
case s = the second line of Eq. ^ is empty. 

If the Donnan potential does not vanish, <f>r> ^ 0, a 
difference in solvation free energy leads to an unequal 
partitioning of cations and anions on the two half spaces 
occupied by solvent A and B. A measure for this unequal 
partitioning is the integrated charge density of the half 
space z < 



(J A ■ = 



dz y^ ag a (z) 
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For KiS <^ 1 the integrated charge density a a is con- 
stant to leading order in KjS, i. e., global quantities de- 
scribing the ion partitioning within the Verwey-Niessen 
model (s = 0) are not influenced by finite interface exten- 
sions smaller than the interfacial Debye length kJ . This 
finding a posteriori justifies the application of the origi- 
nal Verwey-Niessen model to calculate droplet charges in 
Ref. [3. 

The interfacial tension, however, is known to be highly 
sensitive to details of the interfacial structure. In terms 
of the density functional fi (see Eq. (QJ) the interfacial 
tension in excess of the pure, salt-free liquid-liquid inter- 
face is given by A7 = Cl[g+, Q-] - £l[g T ei, ft-ef], where g rci 
is the steplike reference ion number density profile. For 
the excess interfacial tension with respect to the dielec- 
tric interface at z = the reference density is defined by 
g le {(z < 0) := gA and g Te ((z > 0) := gs which leads to 
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(ncosh(Kis) + sinh(Kis))y/gB. 



As the second term on the right-hand side of Eq. © is 
of the order 0(— \/~g~B~) for both gs — > and gs — ► 00 



one finds the following asymptotic behavior of the excess 
interfacial tension: 
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As n and p are experimentally accessible, one can use 
Eq. (J7J) to determine (/>£> or s. The crossover, where the 
low-density asymptotics A7 = 0{— \J~Qb) and the high- 
density asymptotics A7 = 0(±qb) are of the same mag- 
nitude, takes place at the ionic strength 
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For Q B > Qa one finds A^(g B <§C g B ) < and A^(g B > 
p^) > 0, i. c., the excess interfacial tension vanishes 
near the crossover. For qb < Qa, on the other hand, 
^{qb) < for all qb- As the two bulk ion concen- 
trations qa and qb are proportional to each other within 
the present model, one can choose either one calling it the 
ionic strength I. Equation ((SJ) leads to a corresponding 
crossover ionic strength I x . 

Equations and ((HJ) are the main results of the 
present work which will be discussed in the following. 

The results presented so far have been derived from 
the linear Poisson-Boltzmann equation Eq. (J3]) which is 
expected to be reliable if |^>(.z)| <C 1, i.e., \<J>d\ <C 1- How- 
ever, upon solving the non-linear Poisson-Boltzmann 
equation derived from Eqs. ([1]) and ([2|) numerically, we 
found the same asymptotic dependence on the ionic 
strength I, A 7 (7 < 7 X ) = 0(— /?) and Ay(7 > I x ) = 
C(±J), as in Eq. with a crossover at I x > I x where 
the difference I x — I x increases with \<j>r>\- Hence the 
asymptotic scaling of the interfacial tension difference A7 
with the ionic strength I and the existence of a crossover 
I x are robust qualitative features of the linear the- 
ory when compared to the non-linear Poisson-Boltzmann 
theory. Moreover, by numerical fitting one obtains renor- 
malized parameters <j)* D and s* in Eq. ^ such that A7 
calculated within non-linear Poisson-Boltzmann theory 
is reproduced even quantitatively. 

As the asymptotic behavior of the excess interfacial 
tension A7 = 0(±I) for I > I x in Eq. |0 involves 
the parameter s one concludes that the finite size of 
the interfacial region is responsible for this asymptotics. 
This finding is confirmed by published measurements of 
liquid-liquid interfacial tensions |6| and is in fact well- 
known from liquid-gas surface tensions [3, 8]. In con- 
trast, the behavior A7 = 0(-y/T) for K I x in Eq. 
can be attributed to the unequal ion partitioning be- 
cause the prefactor of the asymptotics contains a term 
of electrostatic origin which vanishes if <f>n = 0. The 
latter regime, which gives rise to a negative contribu- 
tion to the interfacial tension, is in contradiction to the 
Onsager- Samaras limiting law C(— 71n(/)) [l0|, which 
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FIG. 1: Magnitude of the excess interfacial tension | A.7I as a 
function of the ionic strength I in water of a water-decaline 
interface for the three salts KSCN (A), KC1 (O), and KI 
(■) according to Ref. 0. The interfacial tension of a salt- 
free water-decaline interface is 50.94 mNm" 1 . The dashed 
lines are power laws ~ 7 whereas the solid line is a power law 
~ I ' both derived from the present model (see Eq. (0). 



contributes positively [ll|. However, according to the 



present model, the image charge interaction is neglected 
in comparision to the electrostatic potential due to the 
unequal ion partitioning, whereas it is the dominating in- 
teraction within the Onsager-Samaras model [lfj . There- 
fore it can be concluded that unequal ion partitioning, 
which is expected to be a general phenomenon for liquid- 
liquid interfaces [12 |, leads to A7 = 0{—\fl) for small J, 
whereas the absence of unequal ion partitioning gives rise 
to A7 = 0(-Iln(I)) 0- The situation of a liquid-gas 
surface with non- vanishing ionic strength in the gas phase 
investigated by Nichols and Pratt [ll| can be considered 
as the borderline between both scenarios such that fea- 
tures of both, the square root and the Onsager-Samaras 
limiting law, can be visible. 

From Eq. ((8]) one infers a high sensitivity of the 
crossover ionic strength I x from the low ionic strength 
regime j(I <C I x ) = 0(—y/l) to the high ionic strength 
regime 7(1 ^> I x ) = 0(±I) upon the model parameters 
(j)D, s, and p, i. e., upon the material parameters of the 
system. This obervation is also borne out by the results 
of Onuki [l2j. Hence, depending on the actual system 
under investigation, I x can be larger or smaller than the 
experimentally available range of ionic strength as will 
be shown in the following. 

Figure Q] displays the magnitude of the excess interfa- 
cial tension | A.-y| of a water-decaline interface as a func- 
tion of the ionic strength I in water for three different 
salts, KSCN (A), KC1 (O), and KI (■), as published 
in Ref. [5|. The interfacial tension of a salt-free water- 
decaline interface is 50.94 mNm -1 . The dashed lines are 
power laws ~ I passing through the largest data points 
for KSCN and KC1, whereas the solid line is a power 



law ~ I 1 / 2 passing through the smallest data point for 
KI. Within the present model one concludes from Fig. Q] 
that the crossover ionic strength I x for KSCN and KC1 
is smaller than 0.01M whereas for KI it is larger than 
0.5M. The prediction A 7 (7 < I x ) < from Eq. © is 
in agreement with the data for KI in Ref. [5(. Finally, 
the excess interfacial tension measured in Ref. [5J is neg- 
ative for KSCN and positive for KC1. Within the present 
model this observation is to be interpreted as follows: 
For the case of KSCN (A 7 (7 > I x ) < 0) one infers 
p > 1 from Eq. (J7J), and consequently in this case solvent 
A is water and solvent B is decaline, because the ionic 
strength in water is larger than in decaline. Assuming 
p > 1 Eq. leads to A-/(I > I x ) ~ -2sl with I = g A 
which, for KSCN, yields s rs 0.53 nm. Hence the discon- 
tinuity of the solvent-induced potentials V a for KSCN is 
located at a distance 0.53 nm on the decaline-side of a 
water-decaline interface. For KC1 (Aj(I ^> I x ) > 0), on 
the other hand, p < 1 due to Eq. (O, i. e., here solvent 
A is decaline and solvent B is water. Assuming p -C 1 
gives rise to A 7 (7 ^> I x ) ~ 2s/ with I = qb which, for 
KC1, leads to s rs 0.23 nm. Thus the discontinuity of 
14, for KC1 is located at a distance 0.23 nm on the wa- 
ter-side of the water-decaline interface. These findings 
suggest a weaker affinity of CI - for the organic decaline 
phase than [SCN]~, which agrees with the structure of 
these anions. Hence the excess interfacial tension data in 
Ref. [5j can be consistently described in terms of Eq. ([7]) 
with respect to the sign and the power law in the ionic 
strength. Moreover, s is, as expected, comparable to the 
size of the ions. However, a more detailed experimental 
check would be highly appreciated. 

To conclude, it has been found within a simple model 
that at small ionic strength / the excess liquid-liquid 
interfacial tension of electrolyte solutions behaves as 
0{— VT) due to an unequal partitioning of ions, whereas 
at large ionic strength it behaves as 0(±I) due to a fi- 
nite interfacial thickness. These asymptotic regimes are 
in agreement with the findings of Nichols and Pratt 11 1 
and Onuki [12]. The crossover strongly depends on the 
components of the system such that all suggested asymp- 
totic regimes can be realized experimentally by choos- 
ing appropriate liquids and electrolytes (see Fig. Q] and 
Ref. [5J ) . The decrease of the liquid-liquid tension at low 
I is expected to be much more pronounced when highly 
charged colloids are considered instead of low-valency 
ions. Quantitative understandig of this increase is of di- 
rect relevance for the stability of Pickering emulsions as 
mentioned at the beginning of this work. On the basis 
of the present simple model for liquid-liquid interfaces 
between electrolytic solutions, which in this work has 
proved to agree with experimental data, it should now 



be possible to study the effect of adding surfactants and 
colloids on the interfacial tension in order to ultimately 
obtain a fully microscopic theory of the (in)stability of 
emulsions [2J, [l2| . 
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